Mon. Not. R. Astron. Soc. 000, 000-000 (2003) Printed 2 February 2008 (MN M?eX style file vl.4) 



The effect of differential refraction on wave propagation in 
rotating pulsar magnetospheres 



D. Fussell and Q. Luo 

School of Physics, The University of Sydney, NSW 2006, Australia 



o 
o 

(N 
oo 



> 
O 

o 
m 
o 

o 

Oh 

6 



Received in original form October, 2003 



ABSTRACT 

Refraction of wave propagation in a corotating pulsar magnetospheric plasma is consid- 
ered as a possible interpretation for observed asymmetric pulse profiles with multiple 
components. The pulsar radio emission produced inside the magnetosphere propagates 
outward through the rotating magnetosphere, subject to refraction by the intervening 
plasma that is spatially inhomogeneous. Both effects of a relativistic distribution of 
the plasma and rotation on wave propagation are considered. It is shown that refrac- 
tion coupled with rotation can produce asymmetric conal structures of the profile. 
The differential refraction due to the rotation can cause the conal structures to skew 
toward the rotation direction and lead to asymmetry in relative intensities between 
the leading and trailing components. Both of these features are potentially observable. 
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1 INTRODUCTION 

Pulsar radio emission is thought to originate in the re- 
gion deep inside the pulsar magnetosphere populated with a 
relativistic electron-positron pair plasma (e.g. Blaskiewicz, 
Cordes, & Wasserman 1991; Melrose 2000). The radio waves 
propagate through the magnetosphere, subject to reflection 
and refraction (e.g. Barnard & Arons 1986; Petrova 2000; 
Fussell, Luo & Melrose 2003, and hereafter FLM) or absorp- 
tion (Blandford & Scharlemann 1976; Lyubarskii & Petrova 
1998; Luo & Melrose 2001; FLM). Since these propagation 
effects can give rise to observable features in the pulse pro- 
file and polarization, study of them can provide insight to 
the radiation processes inside the magnetosphere. 

In the polar cap model, relativistic electron positron 
plasmas are produced in the cascade above the polar 
cap (PC) as a result of rotation-induced particle accel- 
eration (Sturrock 1971; Ruderman & Sutherland 1975; 
Arons & Scharlemann 1979; Daugherty & Harding 1982; 
Zhang & Harding 2000; Hibschman & Arons 2001). Pri- 
mary electrons (positrons) are accelerated to ultra-high en- 
ergies by a rotation-induced parallel electric field, emitting 
high energy photons through curvature radiation or inverse 
Compton scattering. High energy photons decay into elec- 
tron/positron pairs in the strong pulsar magnetic field, form- 
ing an outflowing relativistic pair plasma, referred to as the 
pulsar plasma, which has a very broad distribution in paral- 
lel momenta. Although the mechanism for the radio emission 
is not well understood, it is generally believed that the emis- 



sion is produced from collective plasma processes in the pul- 
sar plasma, through either maser emission or plasma insta- 
bilities (e.g. Melrose 2000). Regardless of the specific emis- 
sion mechanism, the radiation must be produced directly on 
or converted to modes that can escape the magnetosphere 
to reach the observer. In the pulsar plasma, there are two 
escape modes, the X mode with polarization perpendicular 
to the local plane of the magnetic field and wave vector, 
and the LO mode with polarization in this plane (e.g. Ken- 
nett, Melrose & Luo 2000, hereafter KML, and references 
therein). Observations of jumps between two orthogonal po- 
larizations in the position angle of the radiation support the 
hypothesis that the radio emission propagates through the 
magnetosphere in these two orthogonal modes (Stinebring et 
al. 1984; McKinnon & Stinebring 2000). Therefore we only 
consider propagation of these two modes. 

A pulsar plasma is spatially inhomogeneous due to the 
radial and transverse dependence of the density in the open 
field line region. Wave propagation in such an inhomoge- 
neous medium can be affected by refraction and refiection, 
which ultimately determines the emergent pulse profile. Be- 
cause the dipole field, B, decreases in the radial direction as 
B oc where R is the radial distance to the star's cen- 

ter, the density of the plasma that streams along the field 
lines must have the same radial dependence, N oc 
corresponding to the longitudinal characteristic length scale 
for the density variation, L ^ R, considerably larger than 
the star's radius. Variation in the plasma density in the di- 
rection perpendicular to the field lines is mainly due to the 
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electron/positron pair cascade being nonuniform across the 
PC, which gives rise to a much smaller (than R) transverse 
length scale of order the radius of the open field line cone. Al- 
though the inhomogeneities can be on a still much smaller 
scale, possibly as a result of nonstationary acceleration or 
highly localized pair cascades, we only discuss the case that 
both radial and transverse scales are much larger than the 
relevant wavelength and that the refraction is the domi- 
nant effect that changes the ray path. Wave propagation 
can be then treated in the geometric optics approximation 
in which two orthogonal modes propagate independently of 
each other (Barnard & Arons 1986). 

The effect of refraction on wave propagation in the pul- 
sar magnetosphere was discussed in the cold plasma ap- 
proximation by several authors (e.g. Barnard & Arons 1986; 
Petrova 2000; Weltevrede et al. 2003). Numerical study of 
electron/positron pair cascades above the PC shows that 
the distribution of a pulsar plasma has a relativistic spread 
in the plasma rest frame (e.g. Daughcrty & Harding 1982; 
Zhang & Harding 2000; Hibschman & Arons 2001; Arendt & 
Eilek 2002). The relativistic distribution strongly modifies 
the wave properties (KML) and therefore a self-consistent 
model for wave propagation needs to include the relativistic 
effect of the plasma. Past work on refraction of wave prop- 
agation generally ignores rotation and the ray path initially 
in the field line direction would remain two dimensional, 
confined to the plane of the magnetic field lines. The resul- 
tant pulse profile has the same symmetry as the intensity 
distribution at the emission origin (e.g. Petrova 2000). Ob- 
servations often show pulse profiles with asymmetric multi- 
ple components, i.e. there is asymmetry in intensity as well 
as location in pulse longitude of conal components (Lyne 
& Manchester 1988; Gangadhara & Gupta 1998; Gupta & 
Gangadhara 2003) . It is believed that distortion of the pro- 
file is due to aberration (Blaskiewicz, Cordes & Wasserman 
1991; Hirano & Gwinn 2001; Gupta & Gangadhara 2003) 
or absorption (Luo & Melrose 2001; FLM). Wave propaga- 
tion in rotating magnetospheres was recently discussed in 
FLM with emphasis on the cyclotron absorption. Cyclotron 
absorption including rotation can lead to difi^erential absorp- 
tion producing asymmetric pulse profiles. In their discussion, 
FLM assumed the resonance region to be at a substantial 
fraction of the light cylinder (the radius at which the coro- 
tation speed equals the light speed, c) where refraction can 
be ignored. 

The purpose of this paper is to explore the effect of 
refraction on wave propagation in the relativistic pulsar 
plasma that corotates with the star and the implication 
of such effect for interpretation of pulse profiles. We con- 
sider effects of both a relativistic distribution of the plasma 
and rotation. Due to rotation rays that originate from the 
leading and trailing components are subject to different re- 
fraction. It is suggested that such differential refraction can 
significantly distort the pulse profile. Strong refraction oc- 
curs in the region with a large density gradient, which is 
referred to as the refraction region and is located well be- 
low the cyclotron resonance region (FLM). To concentrate 
on the effect of refraction, one assumes the strong magnetic 
field approximation, in which the X mode is not affected 
by the medium and propagates approximately in a straight 
line through the magnetosphere. The LO mode is strongly 
refracted and is discussed here. Following a similar proce- 



dure to that described in FLM, the ray path in the rotat- 
ing medium is evaluated numerically in the geometric optics 
formalism. The emergent pulse profile is obtained at radial 
distances beyond which the refraction is no longer effective. 

In Sec. 2, the wave dispersion of relativistic pulsar 
plasma is summarized with emphasis on the LO mode in 
the strong magnetic field approximation; The formalism of 
ray tracing in the rotating magnetosphere is described in 
Sec. 3. The ray path is obtained by numerically solving the 
ray equations including the relativistic distribution and ro- 
tation. The result is applied to interpretation of asymmetry 
of pulse profiles (Sec. 4). 



2 RELATIVISTIC PULSAR PLASMA 

The wave properties can be obtained based on the relativis- 
tic model of a pulsar plasma (e.g. KML). The distribution 
of the secondary pairs is broad, characterised by a cutoff at 
the lower energy near the (pair production) threshold and a 
rapid increase to a peaJj with the Lorentz factor at around 
7 ~ 10^, followed by a decay and then a cutoff at the pair 
producing photon energy of about 7 ~ 10^ — 10'* (Zhang 
& Harding 2000; Hibschman & Arons 2001; Arendt & Eilek 
2002). It is convenient to use the plasma rest frame, i.e. the 
center-of-momentum frame with the transform velocity de- 
fined by Vs = (7")/(7), where 7 = 1/(1 — v^)-''^^, v is the 
velocity (in c), and (...) is the average over the plasma dis- 
tribution. The distribution in the plasma rest frame has a 
much simpler form, which has an approximately symmet- 
ric peak with a relativistic spread and can be approximated 
by the Jiittner distribution (i.e. a thermal distribution with 
a relativistic temperature). In practical calculations of the 
plasma dispersion it is reasonable to use a much simpler 
distribution such as the bell-type distribution (Melrose & 
Gedalin 1999). This latter distribution is adopted here. 

2.1 Relativistic dielectric tensor 

The response of a relativistic plasma is fully determined by 
the three relativistic plasma dispersion functions (RPDFs) 
W{z) = (^-y - z)-^), R{z) = {-,-\v - z)-'), and S{z) = 
{7 '^{v — z) where z = uj/k\'^c is the parallel (to the 
magnetic field, B) phase velocity (in units of c), fcy is the 
projection of the wave vector along the magnetic field. 

As we are interested in the inner magnetosphere region 
where refraction is important, we adopt the low- frequency 
limit, w <§; fie, where Q.e is the cyclotron frequency, and the 
non-gyrotropic approximation. Only 11^(2;) is required in de- 
termining the dispersion. Specifically, in the strong magnetic 
field approximation the dielectric tensor is reduced to (Mel- 
rose & GedaUn 1999; KML) 

2 

/sTii « /sr22 « 1, K33^l-'^z''W{z), (1) 

where the off-diagonal components are ignored, the 3-ax:is is 
along B, the wave vector is in the 1-3 plane at an angle to 
B, Wp = (e^N leova^Y^"^ is the plasma frequency, and N is 
the total number density of electrons and positrons. 

The properties of Wiz) are discussed in detail in KML. 
One of its main features is that the function has a peak 
near z = w/feyc « 1, which is characterized by the intrinsic 
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relativistic spread (7') > 1 (where 7' — 1/(1 — v'^y^^ is 
the Lorentz factor in the plasma rest frame) and is not par- 
ticularly sensitive to the specific form of distribution. We 
consider a normalized soft-bell distribution in the plasma 
rest frame given by 



F'{u) = 



15 {v'^-v'fivL + v'f 
xH{v'm - v')H{v' + v'„i), 



(2) 



where H{x) = 1 for a; > and H{x) = for a; < 0, m' = ^v' 
and -u'^ < u' < u'^ ^ 7m 7m = 1/(1 - v'l.y^^, v'^ is 
the maximum velocity in the plasma rest frame. Using (2) 
W{z) can be derived analytically. Since F'{u') is a Lorentz 
invariant, substituting?;' = {v—Vs)/{l—Vsv) for (2), the cor- 
responding un-normalized distribution in the pulsar frame, 
F{u) = F'{u'), can be obtained as 



F{u) = 



167m «m (1 - f 

XH{Vm,ax - V)H{V - Vmin), 



(3) 



where Vmax = (Vs + v!^)/{l+v'raVs), Vmin = (Vs - l'm)/(l ' 

v'mVs), Ts = 1/(1 — Us)^^^ is the bulk Lorentz factor. 



2.2 The LO mode 

In the strong magnetic field, the dispersion of the X mode 
is close to that of vacuum waves and since this mode is not 
refracted, it is not considered here. The dispersion of the LO 

mode can be written as (Melrose & Gedalin 1999; Melrose, 
Gcdalin, Kennett, Fletcher 1999) 



z'^ — 1 — tan^ 6kb ' 



(4) 



where 9kb is the propagation angle to the magnetic field, B, 
and the RPDF in the pulsar frame is given by 



W{z) 



VmTg '^{Z^ - 1) [8Ws«77r 



16wJii7m(l - z'^Y 
-(1 -I- v1z^){\ - 3t;i) - (3 - v'i,){z^ -I- vl)\ 
-37m'(l + vsf{z - \f [(1 3«S.)(1 - VsZ) 

-(..-.)(3+.^)]in (;-^-)(;+^-^| 

-i7-'(l - v,)\z + 1)^ [(1 + 3t;^)(l - v,z) 

—^lm{Vs - z)(l - VsZ) [Vmil - VsZ)'^ - {Vs - zf 



X In 



1 - 



(1 - VmVs)[{Vm + Vs - z{l + VmVs)] 



(5) 



where Vm, replacing v'^, is the maximum velocity in the 
plasma rest frame. The average Lorentz factor is given by 



(7) 



i5r. 



W-YmV?, 



2^i^ + ^Vm 



15v„ 



15w, 



+4:VmVs + SVmVs) - 4(1 + 2vm + 5v^ - 5Vs 
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Figure 1. Plots of W(z) for Fs = 2 (dashed) and 5 (solid). 
Although the peak location depends on Fs, for given (7'), the 
function remains approximately the same for z > 1. 



o 2 2 



xvl) In 



/ 1 — I'm \ 
\l+Vm) 



(6) 



In the relativistic limit {Vs ^> 1, 7m ^ 1) the average 
Lorentz factor in the pulsar frame can be approximated by 
{7) « 2rs{7'), where (7') « 57m/16 is the average Lorentz 
factor in the plasma rest frame. 

Since v < 1, W{z) is regular a.t z — 1 and hence, the 
superfacial singularity ai z — 1 in (5) can be avoided by in- 
terpolation using W{z) « W{l) + {l-z)W'{l) with W{1) = 
{7"^(1 - w)"^) = (7(1 -I- vf). Figure 1 shows plots of W{z) 
with Fs — 2 and 5 for 7^ = 2 (corresponding to (7') « 1.2). 
In the relativistic limit, one has W{1) ^ 4(7) 8(7')rs. 
The peak value is Wmax(2) > W{1). 

The two relevant characteristic frequencies can be ob- 
tained from the dispersion relation: the cut-off frequency 
of the Langmuir mode at = ujp{j^'^)^'''^ , the lowest al- 
lowed frequency, which can be derived by setting z —f 00, 
and the crossing frequency of the LO mode uico ~ ujp{'y)^^^, 
at which with decreasing plasma density, the wave evolves 
from the superluminal to the subluminal regime by cross- 
ing the light line {z = 1). The LO mode crosses back to 
the superluminal regime above lOco- In the relativistic limit, 
one has (7"^) w 1/(7) and then the cut-off frequency is 
ujc « 0)^/(7)^/^, which is lower than that for the cold plasma. 
In the mildly relativistic regime one has even a lower cut-off 
frequency ojc = u}p{'y~^}^^'^ < ujp/ {'yY^'^ . 



3 RAY TRACING IN ROTATING 
MAGNETOSPHERES 

Consider pulsar plasma corotating with the star and as- 
sume that the relevant length scale for the inhomogeneities 
is much larger than the wavelength. It is further assumed 
that in the corotating frame the magnetic field is a static 
dipole. The dipole approximation may not be valid near the 
light cylinder and as we consider only propagation in regions 
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well inside the light cylinder where deviation from the dipole 
approximation can be ignored. 

3.1 Hamilton's equations 

In the geometric optics approximation, the ray path is de- 
termined by a set of Hamilton's equations for the spatial 

coordinate x and the wave vector k. To include rotation 
the Hamilton's equations are written in the covariant form 
(Weinberg 1962), 

dx^ _ ODm 

ds dk^ 

dkfi _ ODm 

ds dxi^ 

where the metric is g'^" = = diag(l, —1, —1, —1), x'^ = 
{t, x), fe'' = (w, k), Dm = is the dispersion relation for the 
M mode, which is a Lorentz invariant. The solution to (7) 
gives the ray path, represented by a;^(s) and kfj,{s), where 
s is the parameterization of the ray path. The right-hand 
sides of Eq. (7) can be written as 



(7) 



dkf, 
ODm 
dxu 



4i 

6 



-^ib + 2n) , 



dh 

^ dx. 



n + 



d\nN \ 



^)W'{z)-2zW{z)] 



= [2W{z) + z{z'' - l)W'{z)\ + 2n^ 



6 = 



-^{z'-l)W{z), 



(8) 

(9) 

(10) 

(11) 
(12) 



where b = B/_B is the field line direction, n is the refractive 
index, n = k/fc is the unimodulc of the wave vector, and N 
is the local plasma density. Due to the field line curvature 
and density gradients, one has dk^ /ds =^ 0, which causes 
deviation between the wave vector k and the group velocity 
Vg = dyi/dt. As the ray propagates outward the coefficients, 
fi_3 oc Wp, decreases rapidly, at least as When the 

ray reaches a radial distance, referred to as the refraction 
limiting radius (RLR), beyond which the refraction is no 
longer effective, k and Vg merges into the same direction 
again. 



3.2 The efFect of the relativistic distribution 

The relativistic effect is described by W{z) and W'{z) in ^i, 
^2 and ^3. In the cold plasma approximation, we have 



(z^ - l)W{z) 



n{vs-zr 

1 



(1 - n) 



(13) 



Eq. (7) reduces to the familiar ray equations in the cold 
plasma approximation (Barnard & Arons 1986; Petrova 
2000). Refraction due to density gradients is described by 
^3 OC (z^ - l)W{z). For n « 1, one has (z^ - l)W(z) « 
4(^fcf,r.)/[(6'fef,r,)" + 1]2 < 1/r,. in general, the effect of 
refraction decreases for increasing bulk Lorentz factor, Fg. 
Refraction initially is small at 9kb ^ 1 and rapidly increases 
to reach its maximum at 9kb ~ 1/Fs. 
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Figure 2. Deflection of rays due to the radial gradient. The soUd 
curves (from top to bottom) corresponds to (7') = 5, 10 and 30. 
The refaction of rays strongly depends on the spread of the dis- 
tribution. Compared to the cold plasma case (the dashed curve), 
strong refraction occurs in a much narrower frequency range to- 
wards the cut-off frequency. The bulk Lorentz factor is F^ = 30. 




Figure 3. Deflection of rays due to the transverse gradient. The 
curves from top to bottom correspond to frequencies u) = 10® — 
]^09.76g-r rpj^g ladial distance, R, is in units of the stellaj: radius 
R*. Other parameters are Fs = (7') = 100, Ro = 25, ei = £2 = 3. 



For an intrinsically relativistic plasma, W{z) is strongly 
peaked near z « 1 with Wmax(«) ~ 8(7')F3 (cf. Figure 
1). Near the peak, assuming z ^ Vs one may estimate 
{z'^-l)W{z) ~ iy(l)/8F?{7')2 « 1/F,(7'>, which is smaller 
than that for the cold plasma by a factor of (7'). Inclu- 
sion of the relativistic distribution tends to suppress refrac- 
tion in general and shifts the strong refraction regime to 
u! < Uco, with the strongest refraction found near the cut- 
off frequency, ujc.- Let 6 be the colatitude of k, which is the 
propagation angle with respect to the magnetic axis. A com- 
parison of refraction effects for different (7')'s is shown in 
Figure 2. The ratio of the final (deflected) propagation an- 
gle. Of, to the initial propagation angle, 6i, is plotted as a 
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function of wave frequency ui (in lJc = Wp/ {yV^^)- The den- 
sity is assumed to have a radial gradient oi N ex l/R^ with 
a uniform transverse distribution. It can be seen that the 
effect of refraction becomes negligible for uj/lOc > 1. 



3.3 Transverse distribution of the plasma density 

Strong refraction arises from the effect of the transverse gra- 
dient. The major contributing term to Eq. (9) is from the 
density gradient, dlnN/&x, in the direction perpendicular 
to the field lino, duo to the nonuniform pair cascade above 
the PC. In the current PC model, significant pair production 
occurs near the last open field lines that have the smallest 
radius of curvature. There are few pairs produced on the 
field lines near the polar axis and no pairs on the last close 
field lino on which the accelerating electric field is zero. The 
magnetic dipole field lines in the polar coordinate in the 
corotating frame can be described hy R = sin^ x/sin^ x*) 
where x* ^ Xd is the magnetic colatitude, Xd = {iiR*/cy^^ 
is the half open angle of the PC, and CI = 2ir/P, P is the 
pulse period. The radius of the open field line cone as a 
function of radial distance is R± = XdR^^^, which is the 
characteristic length scale of the transverse gradient. 

For a plasma streaming along the magnetic field lines, 
the transversely spatial distribution of the plasma density 
at a given radial distance is completely determined by the 
density profile at the PC. The density profile is described by 
a function of the magnetic colatitude, X: which is essentially 
the distribution of pair production across the PC. The PC 
is assumed to be circular in the corotating frame with the 
opening half angle Xd- Assuming that pair cascades have an 
azimuthal symmetry relative to the polar axis, the plasma 
density in the corotating frame is peaked at the field line 
with the colatitude on the PC given by Xc ~ £cXd with 
Ec < 1- Specifically, the density is written as (Barnard & 
Arons 1986; Petrova 2000) 



exp 



exp 



n £1 J 



2p 



(14) 



where 'I' — x/XcR^^'^ < I/ec, £1,2 are two parameters that 
characterize the steepness of each side of the peak, is the 
density at the star's surface, the radial distance is in i?*. In 
polar models, the density can be written as = MNqj, 
where Nqj is the Goldreich-Julian density and M is the 
pair multiplicity. Since the maximum pair production oc- 
curs near the last open field lines, one has Xc ~ Xd- Hence, 
it is reasonable to assume that £2 > ei- The Gaussian dis- 
tribution p = 1 in the nonrotating case was discussed by 
Barnard & Arons (1986), also Petrova (2000). We consider 
a more general case including higher density gradients p > 1. 

Figure 3 shows ray paths in the 6-R plane. The rays are 
emitted at -Ro = 25 with the colatitude d = —5°, whore the 
negative angle corresponds to the trailing side of the mag- 
netic axis, and they are refracted by the transverse gradient. 
Notice that refraction is sensitive to frequency: the lower fre- 
quency the stronger refraction. Because of strong deflection 
the first three rays (from top) cross the magnetic axis into 
the leading side. 



3.4 Rotating medium 

Rotation can be included in ray tracing as follows. Assuming 
that the wave is emitted in the field line direction in the 
corotating frame in which the magnetic field is a dipole by 
assumption. The ray path is calculated as a function of time t 
in the observer's frame. Specifically, this involves solving the 
Hamiltonian equations with a sequence of transformations, 
first rotation and then the Lorentz transform, 



T = T2{a)T3{nt)T2{ay 



(15) 
(16) 



where the 3-axis is along the pulsar spin axis (ft), the mag- 
netic pole is assumed to be in the 1-3 plane at t = 0, T2 
and arc the rotation with respect to the 2- and 3-axis, 
respectively, a is the inclination angle (relative to the rota- 
tion axis), Vr = f2 X R is the co-rotation velocity at the 
radial distance, 7?, f2 is the angular velocity of the rotation, 
and I/(Vr) is the Lorentz transform. Eq. (15) represents a 
combined rotation, i.e. an anticlockweise rotation about the 
2-axis by a so that the magnetic pole is aligned with fl, 
then rotation about the 3-axis by fit, followed by a clock- 
weise rotation about the 2-axis by a. In the observer's frame 
we have 



ds 
dk" 
ds 



(17) 
(18) 



Ray trajectories can be obtained by numerically integration 
of Eq. (17) and (18) through the refraction zone, the region 
between the emission radius to the RLR. At the RLR, the 
propagation time, t, and the final wave vector, k/, are deter- 
mined. A rotational transform, determined by t — to, where 
to is the time that a vacuum wave would take to propagate 
through the region, is then applied to k/. One can calculate 
the final wave vectors for all rays and the distribution of a 
bunch of rays can be determined at a particular instance. 



4 APPLICATION TO PULSARS 

The specific mechanism for radio emission is not well un- 
derstood. One possibility is that the radiation is produced 
through linear wave-particle interactions. This is only pos- 
sible for subluminal waves (with phase velocity less than c). 
The LO mode can be subluminal only at frequencies near 
the cross frequency lUco within a small range of propagation 
angle, dkb < Om » {2{i)Y''^u)p/Q.e < 1 (Melrose & Gedalin 
1999; KML). The radio emission can be produced in other 
types of wave and then converted to propagate modes. In 
this case it is possible that the radio emission is in the su- 
perluminal modes with frequency below the cross frequency 
uj < (jOco- If this is the case, refraction may have a significant 
effect on wave propagation. 



4.1 Emission beam 

Since the intensity distribution in the emission region is not 
known, one considers only the case of emission from a par- 
ticular height with a 2-dimensional distribution in intensity 
(e.g. Petrova 2000; ELM). The initial intensity profile in the 
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Figure 4. Asymmetry due to difTercntial refraction caused by 
star's rotation. The solid and dashed curves represent respectively 
the ray paths that orginate from the leading and trailing sides of 
a dipole with the inclination a = 7r/2. The rotation is in the 
direction of positive 6. The parameters are as in Figure 3 but 
with the emission radius Ro = 50; p = 2; (P, P) = (0.1 s, 10~^^); 
u> = 10^-^ s~^; the pair multiplicity is assumed to be M = 10*. 




Figure 5. Differential refraction as in Figure 4 with strong de- 
flection. The asymmetry is easiest to be seen in strongly deflected 
rays. The emission radius is Rq = 20. 



co-rotating frame is assumed to follow the density profile 
and be Gaussian in colatitudes with axisymmetry relative 
to the pole, that is, 



Iq{6q, Ro) = exp 



1 



(19) 



1 /2 

where e — {€i,€2) are defined in (14), ^-o = Xo/x<=Ro < 
l/sc, Xo is the colatitude of the emission point. We have 
Xo ~ 2^0 /3, where do is the corresponding colatitude of k 
at the emission origin. Without refraction the intensity dis- 
tribution (19) would simply lead to a single omission cone. 
Because of refraction the distribution of rays changes as they 
propagate away from the emission region up to the RLR be- 
yond which the refraction is no longer important. The emer- 



gent beam is obtained by calculation of the ray distribution 
at the RLR. 



4.2 Asymmetry in pulse profiles 

Many pulsars have multicomponent pulse profiles, com- 
monly characterized by conal structures (e.g. Lyne & 
Manchester 1988; Rankin 1993) The conal components ap- 
pear to shift toward the rotation direction, which is inter- 
preted as the effect of aberration or retardation (e.g. Blask- 
iewicz, Cordes & Wasserman 1991; Gangadhara & Gupta 
2001; Gupta & Gangadhara 2003). A special emission ge- 
ometry such as nested emission cones is often evoked to ex- 
plain the conal structure of the profile. Although pair pro- 
duction above the PC is expected to be rather nonuniform it 
is not clear how such a regular nested conal structure forms. 
In observations, conal components generally show asymme- 
try in intensities, i.e. some pulsars have a stronger leading 
component, while some others have a stronger trailing com- 
ponent. Our model suggests that diflferential refraction can 
produce such asymmetric pulse profiles without appealing 
to the nested conal structure at the emission origin. 

Figures 4 and 5 show asymmetry in refracted ray paths 
in the observer's frame. One assumes the density profile with 
p = 2 and (ei , S2,£c) = (0.64, 0.6, 0.8). The three pairs (solid 
and dashed) of ray path originate inside the critical field 
lines (x» = Xc) and are initially symmetric on emission in 
the corotating frame, with the intensity distribution given 
by (19). In the co-rotating frame, rays emitted at the same 
colatitude on the opposite sides of the pole are subject to 
different plasma gradients, and due to the differential re- 
fraction their paths are no longer symmetrical. The rays 
originating on the leading side of the pole are focused back 
towards the center earlier on their paths than the ones on 
the trailing side. If one decreases the inclination angle a, the 
rotation effect is reduced. The ray paths in Figures 4 and 5 
appear to focus since they originate inside the critical field 
lines. If the density distribution is strongly peaked near the 
last open field lines and all radiation is produced inside the 
critical field lines, refraction can lead to a narrower pulse 
width. 

The beam intensity profiles in the observer's frame are 
shown in Figure 6 for the plasma density profiles of p = 1 
and p = 2. The intensity profiles are obtained by evalu- 
ating the ray distribution at the RLR. The longitude (in 
degree) is the pulse phase, with the zero pulse phase cor- 
responding to the line of sight in the plane of the mag- 
netic axis and rotation axis. The positive longitude corre- 
sponds to the leading side, which is opposite to the defini- 
tion adoped in Blaskiewicz et al (1991). In their Figure 3, 
the negative phase corresponds to the leading direction and 
the zero pulse phase is where the position angle variation is 
the steepest. We assume oj = 10^ s"\ M = 10*, P = 10"^^, 
Ro = 20 and {£i,£2,£c) = (0.4,0.35,0.8) for p = 1 and 
(ei,e2,ec) = (0.64,0.6,0.8) for p = 2. The conal beam is 
split into two nested cones due to the refraction effect, with 
emission produced inside Xc bending towards the magnetic 
axis forming the inner cone and emission produced outside 
Xc deflected away forming the outer cone. Apart from two 
main nested cones, some more complicated conal structures 
at low-level intensities can also be seen. When the refraction 
is strong, as in the case of steep density gradient (p = 2), 
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Figure 6. Beam intensity profiles. The latitude and longitude (in 
degree) correspond respectively to tlic viewing angle (relative to 
the rotation axis) and the pulse phase. The upper: p = 2, V s ~ 
(7) 16, Q = 60°, P = Is. The lower: p = 1, Ts (7) Rs 100, 
a = 90°, P = 0.1s. The zero pulse phase corresponds to that the 
line of sight is in the plane of the magnetic axis and rotation axis. 
The positive longitude corresponds to the direction of rotation. 



the rays forming the inner cone cross the magnetic axis pro- 
ducing bifurcation of the core. Examples of multiple compo- 
nent profilos resulted from the bifurcation are shown in Fig- 
ure 7. The profiles are obtained from the emergent emission 
beam in the upper subfigure in Figure 6 for throe different 
viewing angles (the angle between the line of sight and the 
rotation axis). The first profile, which is obtained for the 
viewing angle i = a = 60°, has a dominant core structure 
as the result of rays being focused toward the magnetic pole 
axis. All three profiles have nested conal structures. Since 
rays originating on the leading and trialing side are subject 
to different refraction, the distribution of conal components 
skews toward the rotation direction with uneven intensities 
for the leading and trailing components. In the absence of 
refraction, the aberration would advance the emission to- 
ward the leading direction (positive longitude) by a phase 



of only about Rq/Rlc ~ 2.4 degree for iio = 20. The differ- 
ential refraction appears to shift the emission further in the 
leading direction (cf.Figure 6). In our ray tracing model, the 
ray path is not only affected by the aberration but also by 
differential refraction. The latter effect is dominant for the 
emission from low altitudes. 



4.3 Frequency dependence 

As in the cold plasma model, refraction strongly depends 
on frequency and is significant at lower frequencies relative 
to the plasma frequency. The observational implication de- 
pends on emission models in which the wave frequency can 
be related to the plasma parameters. If one assumes that 
the radiation frequency is related to the relativistic plasma 
frequency, say w = riLOc, where > 1 is a model-dependent 
parameter. The frequency dependence is then given by a 
function of lJc, which is mainly determined by the model 
for the density gradient. For the transverse gradient model 
discussed in Sec. 3.3, one has stronger refraction at higher 
frequencies (corresponding to a higher ricOc)- This is because 
the emission region with high cjc must have a narrower den- 
sity profile with a larger gradient and the ray is then subject 
to a stronger refraction (e.g. Petrova 2000). 

Although the hypothesis of emission at the frequency 
fixed relative to the plasma frequency predicts a radius-to- 
frcquency relation, u) ~ R^^^^ , roughly consistent with ob- 
servations (Petrova 2000), there is no widely accepted model 
for the emission mechanism that predicts such a frequency 
relation. This assumption requires that the radiation to be 
produced at a variable range of radial distances, which would 
complicate the frequency to emission radius relation (Wel- 
tevrede et al. 2003). 

4.4 Cyclotron absorption 

A ray that propagates outward eventually reaches the cy- 
clotron resonance region and cyclotron absorption occurs 
when the wave frequency equals the cyclotron frequency in 
the electron rest frame. In general the resonance region is lo- 
cated well above the RLR. Whether the resonance zone is lo- 
cated within the light cylinder depends on the bulk Lorentz 
factor of the plasma as well as the relativistic spread. It was 
shown by several authors that the cyclotron resonance condi- 
tion can be satisfied inside the light cylinder for some pulsars 
(Blandford & Scharlemann 1976; Mikhailovskii et al. 1982; 
Luo & Melrose 2001; FLM). If cyclotron absorption occurs 
it can distort the pulse profile through differential absorp- 
tion (Luo & Melrose 2001; FLM). Since the optical depth 
for the absorption is r oc 6^N, rays that originate on the 
leading and trailing sides propagate along different paths 
with asymmetry in 6 and the density, leading to differential 
absorption. 



5 CONCLUSIONS 

The effect of refraction on wave propagation in rotating mag- 
netospheres is considered including both effects of relativis- 
tic distribution and rotation. Since the relativistic plasma 
dispersion function is peaked at around « w 1 and remains 
small otherwise, inclusion of the intrinsic relativistic effect 
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Pulse phase 

Figure 7. Pulse profiles for three different viewing angles for the 
emission beam in the upper panel of Figure 6. a: j = 60° ; b: 55° ; 
and c: 53°. The intensities in relative scale (with the maximum 
intensity in each case being normalized to 1) are plotted against 
the pulse phase (degree). 



tends to suppress of refraction at frequencies uj > ujc- Refrac- 
tion is sensitive to the plasma density profile in the trans- 
verse direction and causes rays to focus and bifurcate, which 
is qualitatively similar to the result from the nonrotation 
case (e.g. Petrova 2000). The focusing effect tends to pro- 
duce a core component even when the emission region has a 
conal geometry. The bifurcation leads to split of the emission 
cone into two or more nested cones, giving rise to profiles 
with multiple components. 

One of the distinct features of the model discussed here 
is the prediction of asymmetry in the pulse profile. The 
conal components are skewed toward the rotation direction, 
which is qualitatively in agreement with recent observations 
(e.g. Gangadhara & Gupta 2001). The predicted profiles also 
show asymmetry in relative intensities between the lead- 



ing and trailing components, which are common in observed 
pulse profiles. It is suggested that similar asymmetry seen in 
observations can be produced by combination of aberration 
and the differential refraction due to rotation. The latter is 
the dominant cause for the asymmetry if the radiation is 
produced at low altitudes. 

Two other modes that are not discussed here include the 
X mode and low-frequency Alfven mode. In general, pulsar 
plasma is gyrotropic due to that the electron and positron 
distributions are not identical. Although the X mode is not 
purely transverse, with a small longitudinal component of 
along the magnetic field, one expects refraction of the X 
mode to be weaker for the LO mode. Inclusion of the X 
mode would lead to pulse profiles of the two modes that arc 
displaced with one mode dominating any particular part of 
the pulse longitude. The observational implication of such 
a displaced profile needs to be explored. The low-frequency 
Alfven mode can be refracted in the low-frequency regime. 
However, the wave becomes subluminal as it propagates to 
underdense regions and is subject to strong damping. 

Pair production above the PC can be oscillatory over 
the characteristic time scale about Ato = ho/c, where ho 
is the typical length of polar gap acceleration. The pulsar 
plasma formed from the cascade may consist of many out- 
fiowing clouds (e.g. Usov 1987; Asseo & Melikidze 1998). 
For ho ~ R*, we have Ato ~ 10"'' s. Since this time can 
be shorter than the ray propagation time, further work on 
the wave propagation including the time-dependence of the 
medium is needed. 
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